UNSTEADY FLOWS OF FLUIDS WITH PRESSURE 
DEPENDENT VISCOSITY IN UNBOUNDED DOMAINS 
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Abstract. In order to describe behavior of various liquid-like materials at high 
pressures, incompressible fluid models with pressure dependent viscosity seem to 
be a suitable choice. In the context of implicit constitutive relations involving the 
Cauchy stress and the velocity gradient these models are consistent with standard 
procedures of continuum mechanics. Understanding mathematical properties of 
governing equations is connected with various types of idealization, some of them 
lead to studies in unbounded domains. In this paper, we first bring up several 
characteristic features concerning fluids with pressure dependent viscosity. Then 
we study three-dimensional flows of a class of fiuids with the viscosity depending 
on the pressure and the shear rate. By means of higher differentiability methods 
we establish large data existence of a weak solution for the Cauchy problem. This 
seems to be a first result that analyzes flows of considered fluids in unbounded 
domains. Even in the context of purely shear rate dependent fiuids of a power-law 
type the result presented here improves some of earlier works. 



1. Introduction 

Navier-Stokes equations, a widely used model describing flows of incompressible 
Newtonian fluids, are usually expressed in terms of the velocity v = (v^, w^, v^) and 
the pressure p in the form 

(1.1) /9*(vj [Vv]v) = -Vp-f z^* Av, divv = 0, 

where p* and z^* are positive constants representing the density and the viscosity of 
a fluid, and v,, + [Vv]v := f + ^^i 

Within the context of continuum mechanics, the system (11. ip appears as a con- 
sequence of the constraint of incompressibility 

(1.2) divv = 0. 
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and the balance of linear momentum 

(1.3) p*(v,t + [Vv]v) = divT, 

where the Cauchy stress T is related to the symmetric part D of the velocity gradient 
Vv through the constitutive equation 

(1.4) T = -pI + 2z/*D. 

Note that due to (11. 2p p is the mean normal stress, i.e. p = — |trT. Setting 
S := 2i/* D it is obvious that S is the deviatoric part of T. 

Stokes, in his famous paper [IB] that can be considered as one of the corner-stone 
of continuum fluid mechanics, carefully discusses the possibility that the viscosity 
of a fluid may vary with the pressure, see also Hron et al. [21] for details. Barus [4], 
and later on Andrade [1] (see also the book by Bridgman [6]) experimentally showed 
that the viscosity grows with increasing pressure exponentially. Further details and 
references to more recent experimental studies can be found in the book by Szeri 
[18] and in the paper by Malek, Rajagopal [31] that reflects the situation before 
2006. Recent paper£| by Bair and Kottke [3] and by Bair [2] report even drastically 
faster dependence of the viscosity on the pressure. It is worth of noticing that even 
at such higher pressures the variation in the density of most liquids, in comparison 
to the variations in the viscosity, is neghgible, as discussed in Rajagopal [SS] or [21]. 
As a consequence, one can model these liquids as incompressible materials fulfilling 
the constitutive equation 

(1.5) T = -pI + 2z/(p)D. 



There are many experimental works going back to observations made by Trouton 
in 1906 [51] and Schwedoff around 1890 (as discussed in [19]) that confirm the de- 
pendence of the viscosity on the shear rate. A representative list of references can 
be found in the most of books on non-Newtonian fluids, see Schowalter [32], Bird, 
Amstrong, Hassager [5], Huilgol [22] or the survey papers by Malek, Rajagopal, 
Riizicka [32] and Malek, Rajagopal [30] . 

There is less available experimental data that characterize the dependence of the 
viscosity both on the shear rate and on the pressure. This is very likely due to 
completely different experimental set-up used for the measurement of the relation 
between the shear rate and the shear stress on the one hand, and for the measure- 
ment of the viscosity-pressure relationships on the other hand. Nevertheless, we can 
provide references where the incompressible fluids with the pressure and shear rate 
dependent viscosity are chosen in order to model high pressure processes in silos 
and journal bearings, and where some support of experimental observations to the 
considered form of the viscosities is available. 

"'^ Further related experimental studies that concentrates on the dependence of v on the pressure 
and the temperature are presented in Casalini and Bair and Harris and Bair |19] . 
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Schaeffer [41j used the model ([112]) -(OD with 
(1.6) (p, |Dp) = ap|D|-i (a>0) 

in modehng and analyzing flows of certain granular materials. Davies and Li [13] 
and Gwynllyw, Davies and Phillips [TS] (see also [23]) considered a model for com- 
putational simulations of pressure and temperature effects in certain viscoelastic 
liquids, in which the viscosity is of the form 

(1-7) 1^ (p, |D|2) = (rj^ + ^^°^|pp-r ) exp(ap) 

with r = 1.46 and a, tiq, r/oo > 0. 

To summarize there are experimental data confirming a (linear, exponential, or 
even faster) dependence of the viscosity on the pressure, and also on the shear rate. 
Thus, neglecting less significant variations in the density of a liquid at high pressure 
one end-up with the model (11. 2p . (11. 3p and the constitutive equation 

(1.8) T = -j9l + z/(p, |D|2) D. 

In what follows, we shall show that (II. 8p is consistent with standard procedures of 
continuum mechanics if one starts with a general implicit relation between T and 
D. 

1.1. Fluids with pressure dependent viscosities within implicit constitu- 
tive theory. A standard 'derivation' of the constitutive equation for a compressible 
Navier-Stokes fluid starts with the assumption that the Cauchy stress T in a fluid 
depends on the density and the velocity gradient Vv, i.e., 

(1.9) T = H(^?,Vv). 

The requirement that T is invariant with respect to any coordinate transformation 
leads to the conclusion (see for example Truesdell [52j or Serrin [43j) that H depends 
on the velocity gradient through its symmetric part D(v) and H is an isotropic 
second order tensor. This means that 

Uig, QDQ^) = QH(£), D)Q^ VQ G Q , 

where Q denotes the orthogonal group. Using the representation theorem for such 
isotropic second order tensors one obtains that 

T = ai{g, Id, //d, ///d)I + «2(£>, /d, //d, ///d)D + a^ig, In, IId, ///d)D2, 
where 

Jd = trD, JJd = trD^ 1 1 Id = trD^ 

If one requires that the stress depend linearly on D, then one obtains the constitutive 
equation for the classical compressible Navier-Stokes fluid, namely 

(1.10) T = -p{g)I + 2fi{g)D + A(^»)(trD)I . 

This straightforward procedure (that derives the constitutive equation for an com- 
pressible Navier-Stokes fluid starting from the assumption (II. 9p and followed by 
some standard requirements) however loses its clarity if one aims to arrive at an 
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incompressible Navier-Stokes fluid. We recall that the classical approach that is em- 
ployed in most continuum mechanics textbooks to enforce internal constraints (such 
as the constraint of incompressibility (11. 2p ) is to require that the constraints do no 
work. Then the possibility that the viscosity depends on the pressure is excluded, 
and one cannot deduce the constitutive equation like (II. 5p or (II. 8p . We refer to 
Rajagopal and Srinivasa [38j for more details and for outlining how to overcome this 
drawback (within the framework of explicit constitutive equations of the form (II. 9p 
in purely mechanical context) by relaxing the requirement that the constraints do 
no work). 

It is thus interesting to observe that the procedure starting with the assumption 
(II. 9p and ending with the compressible Navier-Stokes fluids (ll.lOp is also applicable 
to incompressible Navier-Stokes fluids provided that we start with a general implicit 
constitutive equatioiJl between T and D of the form 

(1.11) G(T,D) = 0. 

Analogously as above, if we require the function G to be isotropic, then G has to 
satisfy the restriction 

G(^,QTQ^,QDQ^) = QG(^,T,D)Q^ VQ G Q . 

The representation theorem relevant to such isotropic tensors immediately implies 
(see Spencer [^) that 

ttol + aiT + aaD + ogT^ + 040^ + a5(DT + TD) + aJT'^B + DT^) 

(1 12) 

+a7(TD2 + D^T) + a8(T^D^ + D^T^) = , 
where a^, i = 0, . . . 8 depend on the invariants 

trT, trD, trT^ trD^ trT^ trD^ tT{TB),tT{T'^B),tT{'D^T),tT{T^B^). 
We note that if 

tti = 1, 

and 

^2 = -2K-^trT, IDH (^>0), 

we obtain (II. Sp . 

^Implicit constitutive theory, as introduced in [2S] and [33 (see also ^S] for further development 
and [26) for a reflection of the implicit constitutive theory on mathematical analysis of non-linear 
PDE's), provides the framework that is sufficiently robust to capture complicated response of 
materials. In this approach the quantities such as stress and strain in solid-like models and stress 
and the velocity gradient in fluid-like models share an equivalent role (on contrary to the explicit 
constitutive theory that usually prefers strain to stress since stress is a function of the strain or 
the velocity gradient). In addition, the implicit constitutive theory can eliminate some internal 
variable theories that are frequently connected with less clear physical meaning and with difficulties 
to identify appropriate boundary conditions. 
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Also, requiring that fll.l2p is linear both in D and both in T, one obtains 

70 (trT) I + 71 (D ■ T) I + 72T + 73 (trT) D + 74(TD + DT) = , 

where 7^, i = 0,1,2,3,4, are constants. Requiring that the constraint of incom- 
pressibility holds, i.e. trD = 0, the last identity implies that (oj, i = 1,2,3, are 
constants) 

ai(T - -tiT I) + a2(-trT)D + a3(TD + DT - 2T ■ D) = . 

3 3 

Setting a3 = we obtain a model where the viscosity depends on p linearly (similarly 
as in the Schaeffer model (11.61) ). 

1.2. On flows in unbounded containers and relevant results established 
earlier. Any real flow takes place in a bounded domain. Nevertheless, in order to 
understand the mathematical properties of the equations that govern fluid motions 
and to observe how solutions relevant to a material A can differ from solutions 
of the system of equations describing flow of another material B one frequently 
comes across unbounded domains. Flows between two infinite parallel plates (plane 
Couette-Poiseuille flows), flows in a tube, flows glowing down an inclined plane, flows 
due to suddenly accelerated plate or flows due to an oscillating plane are examples 
of motions that are considered in unbounded domains. Despite their unboundedness 
the shapes of these domains are very similar to containers where properties of flows 
are experimentally measured. Then the analytical solutions can be compared with 
experimental data and in most cases the coincidence is very good. Several studies 
that aim to obtain some explicit solutions are in place. 

Hron et. al. [2T] investigate the possibility to find explicit solutions for flows 
between two infinite parallel plates with no-slip boundary conditions and for the 
viscosities of two types: (i) z/(p) = exp(ap) and (ii) z/(p, |Dp) = aplDl*""^ for 
r G (1,2). (See also relevant study [47J where however some imprecise statements 
are made, and [50] focused on some stability issues.) Vasudevaiah and Rajagopal [53] 
considered the fully developed flow in a pipe dealing with a fluid that has a viscosity 
that depends on the pressure and shear rate and were able to obtain explicit exact 
solutions for the problem. Kannan and Rajagopal analyzed effects of gravity on 
flows between rotating parallel plates. For various viscosity-pressure relationships 
they observed, among others, that a boundary layer can be adjacent to just one of 
the plate. Massoudi and Phuoc [33] considered flows between parallel plates due to 
oscillatory pressure gradient. Rajagopal [37] finds special solutions for flows down 
an inclined plane. Srinivasan and Rajagopal in [15] study flow of fluids with pressure 
dependent viscosity due to a suddenly accelerated plate and due to an oscillating 
plate. In all these flows the velocity field, and consequently also the structure of the 
vorticity and the shear stresses at the walls are significantly different from those for 
the classical Navier-Stokes fluid. 

In mathematical analysis, in order to separate the problems connected with gov- 
erning equations from those that are due to the presence of the boundary, it is 
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preferable to look for solution either in the whole spac^ or in spatially periodic set- 
ting; the former corresponds to the Cauchy (or initial-value) problem. The Cauchy 
problem for flows of a one class of fluids with pressure and shear rate dependent vis- 
cosity is the topic investigated in this paper. Regarding the approach we carry on 
the regularity method employed by Malek, Necas and Rajagopal |27j for S{p, D(v)) 
and spatially periodic problem, and on Pokorny work [34j dealing with the Cauchy 
problem for S(D(v)), being independent of p. However, we have to strengthen sev- 
eral steps in Pokorny approach in order to make the whole procedure applicable 
to fluid with pressure and shear rate dependent viscosity. As a consequence, some 
results presented here seem to be new even for fluids with shear rate dependent 
viscosity. 

Before going to details we survey results that investigate mathematical properties 
of flows of such fluids in general, but compact domains. 

To date there have been few mathematically rigorous studies concerning fluids 
with pressure dependent viscosity. To our knowledge, there is no global existence 
theory that is in place for both steady and unsteady flows of fluids whose viscosity 
depends purely on the pressure. Previous studies by Renardy [10] , Gazzola [IB] and 
Gazzola and Secchi [T7] either addressed existence of solutions that are short-in-time 
and for small data or assumed structures for the viscosity that are contradicted by 
experiments. Recently, there has been some resurgence of interest in studying the 
flows of fluids with pressure dependent viscosities. Malek, Necas and Rajagopal [28] . 
Hron, Malek, Necas and Rajagopal [20] and Franta, Malek and Rajagopal [15] have 
established existence results concerning the flows of fluids whose viscosity depends 
on both the pressure and symmetric part of the velocity gradient in an suitable 
manner. Franta et al. [12] established the existence of weak solutions for the steady 
flows of fluids whose viscosity depends on both pressure and the symmetric part of 
the velocity gradient, that satisfy Dirichlet boundary conditions. This result has 
been extended to lower values of power-law exponent in [8]. Earlier, Malek et al. 
[27] and Hron et al. [20] established global-in-time existence for unsteady flows of 
such fluids under spatially periodic boundary conditions. The extension of these 
results to flows in bounded domains subject to the Navier's slip are due to Bulicek, 
Malek and Rajagopal [7] and to fully thermodynamical setting in Bulicek, Malek 
and Rajagopal [iOj. There are models where the viscosity goes to oo as p — oo and 
the existence of weak solutions can be established, [9]. The relation of models and 
results to implicit constitutive theory is discussed in a survey paper [2B] . 

2. Setting of the Problem and the Main Result 

We are interested in studying the following Cauchy problem: to flnd (v,p) : 
[0, T] X — )■ X M solving (in a weak sense) 

(2.1) vt + div(v® v) -div(S(p,D(v))) = -Vp, in [0,T]xR^ 

(2.2) divv = 0, in [0,r]xM3 



Leray's fundamental paper [24] on long-time and large-data existence theory for the Navier- 
Stokes equations concerns the Cauchy problem. 
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(2.3) v(0, ■) = vq with div vq = 0, 

where T is the length of the time interval of interest, Vq is given divergenceless vector 
field (initial velocity) and S(p, D(v)) stands for the constitutively determined part 
of the Cauchy stress T being of the form T = —pi + S(p, D(v)). Typically, we 
consider the case S(p, D(v)) = z/(p, |D(v)p)D(v), allowing thus the possibility that 
the viscosity depends on the pressure p and the quantity |D(v)p, that at simple 
shear flows simplifies to the shear rate. 

As regards the structure, we assume that for any 1 < r < 2 there exists Ci, C2 > 
such that for all p G M and all symmetric matrices B, D G W^^'^ 



fAll 



Ci(l + |D|2)^ |B|2 < 



gS(j9,D) 



(B®B) < ^2(1 + IDP)— IBI 



and there is a 70 > such that for all G M and all D G W^^'^ 
aS(p,D) 



(A2) 



dp 



<7o(l + |D|^)— <7o. 



As a consequence of (]Aip - (]A2p we get 

Lemma 2.1. Let assumption (Al) hold. Then for all p eM. and D G 



b3x3 
sym 



(2.4) 

and 
(2.5) 



S(p,D) - D > Ci|D|2(l + |D| 



> — min(|D|MDr') 



S(p,D) 



< CIDKI 



|D| 



< C2min(|D|, |D 



Proof. See Lemma L19 p. 198 in [29j. 

The aim of this paper is to establish the following existence result. 



□ 



Theorem 2.2. Let vq G L\^ and S fulfil the assumptions (]A1|) -( |A2|) with 70 < 



■ Assume that r G (|,2) then there is a couple (v,p) such that 



(2.6) vGL°°(0,r;r 



n L"(0, T; D^'^iR^) H W^'^ (0, T; W' 



■2,2/Tn)3N3 



with (j){s) 



and 



(2.7) pi G L\0, T- L«(M3)), for all q G (1, -), 

o 

(2.8) p2e L'{0,T;L'{W^)), for all s e [2,r'], 

where p is given as p = pi + p2, and that solves for all -0 G L^r-e (0, T; VF^'^(]R^)) 



(2.9) 



[ (vt,V)-(v®v,VV') + (S(p,D(v)),D(V))dt= [ {p, div 1/^) dt. 
Jo Jo 
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Note that all spaces used in the previous theorem are introduced in the next sec- 
tion. 

Nowadays, there are several results concerning the model in investigation. All of 
them are treating flows on bounded domains. The novelty in our study is that we 
consider unbounded domain (especially M^). Regarding the approach we carry on 
the regularity method employed by Malek, Necas and Rajagopal [27] for S{p, D(v)) 
and spatially periodic problem, and on Pokorny work [34j dealing with the Cauchy 
problem for S(D(v)), being independent of p. 

3. Basic definitions and auxiliary lemmas 

In this section we will fix the notations, state the basic definitions and recall 
some known and useful tools. In order to distinguish between scalar-, vector- and 
tensor-valued function we denote by bold symbol the vector-valued function, i.e. 
V := {v^, v^, v^), and capital bold symbol for tensor valued function, i.e., (S)jj := Sij 
for i,j = 1,2,3. Let (/) : M+ — )■ M+ be an increasing continuous convex function that 
vanishes at zero. For arbitrary open C M'^ we denote by the symbol L^{Q) the 
Banach space 

L'^(n) := ^p)"'"^^ with = inf |a > 0; J ) dx < 1 

Note that if (f){s) = with some r e [l,oo), we write := L<^(fi) and this 

definition corresponds to the standard one for Lebesgue spaces. The space L°°{Q) is 
defined by usual way. The Sobolev spaces W'''P{Q) are defined through usual way. 
In addition, for arbitrary A; G Nq and arbitrary </> satisfying the condition mentioned 
above we define the space D^''^{Q) as 

D'''^{Q) := ^'(^])"■"^'■^ with \\v\\Dk.4. := \\V''v\\l^. 

Similarly as above, the symbol D^''^{Q) with r G [l,oo) denotes the space D^''^{Q) 
with (/){s) := s^. Since we work with solenoidal functions, we also denote 

W^r^in) := {vGH/'='W;div^ = 0}, 

where we used the notation X"^ := X x ■ — x X ^. As usually, for a Banach space 

m— times 

X, the symbols U{0,T; X) and C{0,T;X) stand for the standard Bochner spaces. 
Finally, if a G L''(M^) and b & U (M^), we define (a, 6) := Jj^g a{x)b{x) dx, and 
in addition for general Banach space X we also set (a, 6) := {a,b)x,x* whenever 
a & X and b G X* and whenever the meaning of the duality is clear. We should 
note that in the rest of the paper the functions (p, : M+ — t- M+ are defined as 
(j){s) := s^(l + s^)^ and ^{s) := min(s^, s^). 

In addition we recall all auxiliary lemmas needed in the following text. First, 
since we will need to deal with Orlicz spaces, we will frequently use the embedding 
theorem. 
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Lemma 3.1. There there exists a constant C > such that 

(3.1) \\v\\l^ <C\\Vv\\l^ 
for allv e D^''I'(R^). 

Proof. The proof for more general (f) can be found in [121 Theorem 1, p. 432]. □ 

The next lemma is in fact the Korn inequality in some Orlicz spaces. 

Lemma 3.2 (Korn inequality). Let r G (|,2). There exists a constant C > such 
that for all v e D^''t'{R^f there holds 

(3.2) \\Vv\\l,<C\\ D(v)|U,. 

Proof. Since (p is convex functions which satisfies A2 condition and its comple- 
mentarjQ function as well, the inequality (13. 2p is a consequence of [T¥l Theorem 
6.4]. □ 

As a simple consequence of (13.11) and (13. 2p one can deduce the following 

Corollary 3.3. Assume that q G [2, -y]. There is a constant C such that for arbi- 
trary V satisfying 

I|v||l°°(0,T;L2(k3)3) + ||D(v)||^r(o,T;L'/'(R3)3x3) < K 

there holds 

(3.3) ||v||l9(0,T;L'J(R3)3) < CK. 

Proof. First, using (13. ip and (13. 2p . we conclude that ||v||^r(o,r;L'/'(R3)3) < CK. Hence, 
we can estimate (with = ||f ||l'^) 

f \-v\'idxdt= f f |v|'?dx+ f \-v\'idxdt 



< I llv||2^' ( / Ivr dx 

v|<V,;. 



6r — 69+2rg 
I2 



2-2 

4 



3r 



3-r 3r(9-2) 
3r 5r — 6 



+ l|v||,''"-'^ I / |v|^ dx 1 dt 

v|>v^ 



fT 6^ 3(q-2) 6r-69 + 2r9 3r(g-2) 

Jo 

provided that 5 < -y. □ 
Next, we use the convention that the solution to the following Laplace equation 
(3.4) Au = / in M^ 



^The definition of complementary function and also of A2 condition can be found eg. in [T4] . 
Note that our cj) satisfies all the assumptions of the Theorem 6.4 in [T3] because 1 < r < 00. 
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is always given by 

(3.5) u{x) [ dy, 

4vr Jus \x - y\ 

whenever the integral on the right-hand side is meaningful. Using Calderon-Zygmund 
singular integral operator theory, we can conclude that for all g G (1, oo) 

(3.6) \\VM,<K,\\f\U, 

(3.7) \\Vu\\, < K,Mg, with/ = divv, 

(3.8) < A^'^IISIIg, with / = div div S. 

Moreover, it is easy to show that we can set K2 = 1. In addition, the solution 
of f l3.4p given by fl3.5p is unique in the class of function that solves (13. 4p and that 
vanish at infinity. Having theory for Laplace equation, we introduce the so-called 
Helmholtz decomposition 

where for a given v G iy^'^(M'^)^, solves the problem 

(3.9) Ag"^ = div v, g'^ix) — )■ as \x\ — > 00, 

and Vdiv := V — Vg^. The following estimates then easily follows from f l3.6p -f lX7|) 

(3.10) ||VVll,< /^Jldiv v||g, ||Vvrf,,||g < K,||Vv||g, 

(3.11) IIVS'^lls < ^s||v||s, \\Vdiv\\s < Ks\\v\\s, 

whenever the right hand sides make sense. 

Next, we define an orthonormal basis which is essential to perform a Fadeo- 
Galerkin scheme. Recall that for bounded domain such a basis can be formed 
by the eigenvectors of the operator —A with suitable boundary conditions. Unfor- 
tunately, for this procedure fails since the Laplacian has no eigenvalue. Using 
the fact that P(M^) is separable we have: 

Lemma 3.4. There is an orthonormal basis B = {o;^, u;^, ■ ■ ■ } of W^^{M.^) such 
that each uj^ belongs to V(R^). 

Proof. For the convenience of the reader, we outline the proof. Since Pdiv(IR^) is 
dense in W^('^, and r'div(IR^) is separable, there is a countable subset C = {(p^, 0^, ■ ■ ■ } 
of V(R^) dense in W^^;. From C we can construct an orthonormal basis B = 
{a;^, o;^, ■ ■ ■ } of the Hilbert-space W^'^ following the Gram-Schmidt procedure. □ 

Finally, we review the classical Caratheodory theory for ordinary differential equa- 
tions 

(3.12) ^(t)=A{t,m) a-e., ^(0) = x, 

with non-smooth right-hand side ^ : / x — or equivalently the corresponding 
integral equation 

(3.13) C{t) =x + [ A{s,^{s))ds. 

Jo 
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Definition 3.5. Let A : I x M.'^ ^ M.^ be a function, where I G is an interval. 
We say that A satisfies the Caratheodory conditions (CC) if the following holds: 



• A{t, x) is Lebesgue measurable in t for all fixed x G M^, 

• A{t, x) is continuous in x for almost all t ^ I, and 

• sup |^(t,a;)| < (3{t) a.e. for some positive function (3 G L]^^{I). 

Note that the first two Caratheodory conditions ensure Lebesgue measurabihty 
of the composition s — )■ A{s,f{s)) for all / G while the third condition is 

crucial in the existence proof. 

Theorem 3.6 (Existence theorem for ordinary differential equations). Let I be some 
subinterval ofM. and assume that A satisfies (CC ) on I xM.^ . Let x G M.^, then there 
exists an absolutely continuous solution ^ = {^i, ■ ■ ■ , Cd) to the ordinary differential 
equation f l3.12p . 



4. Equation for the pressure 

This section is devoted to the solvability and the properties of the solution to the 
following equation 

(4.1) - Ap = divdiv(v®v-S(p,D(v))). 
The basic properties are then established in the following 

Proposition 4.1. LetS satisfy flAl|) - flA2|) with'jo < 1 andwithr G (1,2). Assume 
that there are 2 < gi < cxd such that v G L''(]R'^)^ for all q G [2, qi], and assume that 
V G D^'''^(R^)^. Then there exists unique p solving (14. ip such that p = pi+P2 where 

(4.2) Ibillf <C(g)||v||; forallqEi2,q,], 

(4.3) \\p2\\s < C{s)\\B{v)\\l, for all s G [2,r']. 




In addition, if ^ E W'^''^{M?Y then the following estimate holds 
(4.4) 



\\Vph < ( C||Vv||2|| Vv||3 + cJ [ (1 + |D(v)|2)'-2|D(Vv 

1 - 7o V V^M3 



Moreover, if there is a sequence {v^jj^-^ satisfying 

(4.5) v" ^ V weakly m L\Q, T; L'^iR^) for all q G [2, q^] 

(4.6) v" ^ V a.e. m {0,T) x 

(4.7) Vv"^Vv a.e. m {0,T) xR^, 

(4.8) II v"!!/)!,,;, < C uniformly w.r.t. n, 
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then there is a (not relabeled) subsequence solving (14.11) such that (p" 
Pi + P2) 

(4.9) pI^Pi weakly in Li{0,T;L2{R^)) for all q e {2,qi], 

(4.10) P2^P2 weakly m L'{0,T; L'{R^)) for all s G [2,r'], 

(4.11) a.e. (0,T) X 

and the couple {p, v) solves (14. ip . 

Although the proof follows almost step by step the procedure developed in 
(see also [7]), we prove it here in all details because in [27117] the similar results were 
proven only in bounded domain. 



Proof of Proposition \4.1\ We prove the existence of p solving (14. ip only for v G 
D(R^)^. The existence result for general v then follows from (l4.2p -f H?TT|l . Hence, 
let V G P(M^)'^ be arbitrary. Define p^ = and we find as the solution of 

(4.12) - Ap" = divdiv (v ® V - S(p"-\ D(v))) . 

Next, we show that is Cauchy sequence in L^(]R'^) and therefore converges strongly 
in L^(M^) to some p that has to solve (I4.12p . To do it, we first note that due to 
flAT]l and G L'^(R^) for all n. Next, since 

-A(p" - p"^) = div div (S(p"-\ D(v)) - S(p'"-\ D(v))) 

the estimate (13. 8p implies that 

1 IA2I 1 

lb"-p'"||2< ||S(p"-\D(v))-S(p™-\D(v))||2 < 7o\\p"-' - p'''-%- 

Since 70 < 1 we end by using the Banach fixed point theorem. 
Next, for a given (p, v) solving (14. ip we define pi, p2 through 

(4.13) -Api = divdiv(v ® v), 

(4.14) Aj92 = divdivS(p,D(v)). 

Therefore we see that pi+p2 = p and (13. 8p implies (after using (12. Sp . the assumption 
on V and the fact that 1 < r < 2) (I4.2p -( HI3|) . Next, using (13. 7p . we immediately 
deduce that 

||Vp||2< |||v||Vv|||2+||VS(p,D(v))||2 

1 

< C||Vv||2||Vv||3 + C2 (^j jl + |D(v)n"t^|D(Vv)|2 di^ ' + 7o||Vp||2 

where for the second inequality, we used standard interpolation and Sobolev embed- 
ding and also flXT|) - fR2D (see [27] for details). Thus, (jOD easily follows. 

Finally, assume that we have sequence v" satisfying (I4.5p - (l4.8p . Hence, we can 
find the corresponding sequence p^ solving (14. ip and satisfying (l4.9p -f H?TU|) . where 



^We use solvability of (|4.1[) at each time t. 
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Pi,P2 are defineci through (KT^ . f imj) . First, we show a.e. convergence of p". We 
find g"^ solving for some fixed s, such that 1 < s < y 

A^7'^ = iK-Pir~2(K-Pi). 

Without loss of generality assume that qi < 3. Therefore we get s < 2. Hence, using 
dSSl) and (gi]) we get that 

qi 



r 2 s 



(4.15) ^ IIV^^^II:^ dt<C for all a e l^- 

Then, we find a sequence of nonnegative functions fk such that fk{x) = 1 in B{0, k), 
fk{x) = in R3\5(0,2A;), and that satisfies |V/| < Ck'^ and |VV| < Ck-^. Next, 
multiplying the equation for — Pi, i.e., 

-A(p" -pi) = divdiv(v" ® v" - v ® v) 

by integrating over (0, T) x and using integration by parts, we find that 

'"^(Pi - Pi, A((7"/fc)) dt=- /"^(v" ® v" - V ® V, V2(^?"/fc)) (it. 



Using, (14. 5 p and (14. 6 p and the fact that we integrate over a compact domain, we see 
that the integral on the right hand side tends to zero as n tends to infinity, provided 
that (14.151) is satisfied for some a > (y)'. Since s < y this is always true. Next, 
using the identity A{g'^fk) = fk^g^ + g^^fk + 2V5f" ■ V/^ and the definition of 
we deduce that 

limsup / / — pi\'^ fk dx dt 

n-5>oo Jo Jr'a 

<limsupc/" [ \p'l-pi\{\Vg''\k'^ + \g''\k-^) dxdt. 

n->oo Jo J B{Q,2k) 

Setting a := 2(s-i) deduce from the fact that gi < 3 that a > 3 and since s < y 

we also obtain that a' < y. Consequently, using flTOj) we obtain ||V^"||^ + 
1 1 9'" I loo — ^- Finally, the Holder inequality, (14. 9 p and the fact that gi < 3 imply 
that 

[ I \pl-pi\{\Vg''\k-^ + \g''\k-^) dxdt<C{k'^ + k-^)k-a Q 

Jo JB(0,2fc) ''^'^ 

and using this it is easy to deduce that — )■ pi a.e. in (0,T) x M^. 

Next, we apply the similar procedure also for jOj- First we find for some fixed 
s G (2, r') for which s < 3, the function g"^"" solving 

^g'^''=\P2-P2r\P2-P2)fk~\ 

where fk is defined above. It is consequence of (13. 6 p and (l4.1Up that 
(4.16) / ||V'^™"||:5<C for all aG^ ^ 







1' s-1 



®The Bochner measurability over time follows from the fact that the mapping that assigns to 
fixed V some p solving (|4.1I) is continuous. We refer to [27j for details. 
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Next, multiplying the equation for — by fug"^^-, integrating over (0,T) x M? 
and using the similar procedure as above, we obtain 

f I \P2-P7\'fk dx dt < r(Sb",D(v")) - S(p™,D(v-)), V^l^?'""/,)) dt 
Jo Jr^ Jo 



T 



+ C I I \p^-p^\{\Vg"'''\k-^ + \g"'''\k-^)dxdt=:I^'' + I^'' 



'0 J B {0,2k) 

To estimate we use flAip - flA2p (see |27] for details) to observe 

ir = Hsip^, D(v")) - Sb^ + P2, D(v)), V^(^/""/fc)) dt 
Jo 

+ [ (S(p^+p2,D(v))-S(p,D(v)),V^(<?™"/,))dt 

+ / (S(p,D(v))-S(pr+P2,D(v)),V2(r"/fc))c^t 

+ / (S(pr+P2,D(v))-S(p'",D(v™)),V2(^7"^"/,))dt 

o(m,n) + 70 / i\P2 ' P2\ + \P2 ' P2l'^\g"'Vk)) dt, 
Jo 



< 



where 



lim o{n, m) = 0. 

n,m—^oo 

Thus, using fl3.6l) and the definition of g^^, we finally derive the estimate 

m-P2)fkrs dt < o{m,n) + 70 T \m-p2)fk\\i + uP2-P2)fkrs dt 

Jo 

Hence, applying for /™" the same procedure as above, taking limit w.r.t. n and 
m, using weak lower semicontinuity of norm and the fact that 70 < |, we finally 
conclude that P2 — )■ P2 a-e. in (0,T) x which consequently implies f l4.1ip . □ 

5. The 5 - approximation 

Throughout this section 5 > is fixed real number. In order to take easily the 
limit from the Galerkin (finite-dimensional) approximation to a continuous (infinite- 
dimensional) approximation we have to define carefully the approximation of (12 .1 p - 
(12. 3p . We introduce the following 5 - approximation: 

(5.1) Vt + div(v®v)-div(S(p,D(v))) + 5A2v = -Vj9 in [0,T]xM^ 



(5.2) divv = in [0,T]x 

(5.3) v(0,-) = (vo)5 in Ml 

Here the symbol stands for usual mollification by convolution. The main exis- 
tence result of this section is the following 



UNSTEADY FLOWS OF FLUIDS WITH PRESSURE DEPENDENT... 15 

Proposition 5.1. Assume that vq G L'^^^{M.^) and 1 < r < 2. Then there exist 
V := v-^ G W^'\0,T; L'^{M.^f)nL'^{0,T;W^{^{M.y), 

and 

satisfying f l5.ip -( l5l2|) weakly, i.e., 

[ (vi, V') - (v ® V, VV') + (S(p, D(v)), D(V)) dt 
Jo 

+ S [ (VV,VV)t^^= / {p,div^i})dt y e L\0,T;W^'\Ry). 
Jo Jo 

Proof of Proposition \5.1[ We proceed via Faedo-Galerkin scheme. We look for = 

N 

v^(t, x) := (t) uj''{x) that is the solution of Faedo-Galerkin approximation of 

k=l 

(l53]) -f l5:2|) .i.e.. we look for a vector c^{t) := (cf (t), c^(t), ■ ■ ■ G that 

solves the following system of ordinary differential equations 

^ ^'^''\u;>'] - (vl®v^,Va;'=) + (S(P(v^), D(v^)), D(a;'=)) 

(5.5) 



(5.4) 



V dt 

+ 5(VV^,V2a;^-) =0, fc = l,---,iV, 
I v^(0) = K^),, 

AT 

where (vi^)^ := ((^0)5, a;''^)iy2.2 a;*^ and 'P(v^) := p^, where is given by (14. ip 

fc=i 

with v^. The initial value problem (15. 5 p can be written in the form 
(5.6) c^ = ^(c^) with c^(0)=C^, 

where J-" : — > is a continuous function. Hence we obtain an ordinary differ- 
ential equation for which the Caratheodory theory can be applied. In conclusion, 
for every > 1 there exists an approximate solution on the short time interval 
[0,r] that can be however extended onto the whole interval of interest (0,T) by 
using the estimates proved below. Our aim is to pass to the limit as — )■ 00 (for 
fixed 6). To do it, we shall first establish a priori estimates that are independent of 
A^. 

Next , multiplying the k-th equation in (I5.5p i by (t), using the fact that con- 
vective term vanishes after integration by parts (v^ has compact support), using 
the estimate (12. 4 p and integrating the result over times interval (0,f), we deduce 
that 



i^dw^'^ml + Cr f I 0(|D(v^)|)rfx.it+ f 6\\V\^\\l<\\\{^^ 
Jo Jm? Jo 



N\ \\l 

J<5||2- 
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Now, usin^ Proposition 14. 1^ fl5.7p and the fact that vq G L^(R'^)^, we deduce the 
inequahtju 



sup ||v^(t)||2 + 2Ci r I 0(|D(v^)|) rfxrft 

E(0,T) Jo ilR3 



(5.8) 



It is worth of noticing that the constant appearing on the right-hand side of (15. Sp 
is 5-independent. Next, we deduce the similar estimate also for time derivative v^. 
Hence, multiplying the k-th. equation in (l5.5P i by (c^(t))f and using integration by 
parts (v^ has compact support), we deduce that 

(5.9) l|vf(t)||^ + ^|||V^v^(t)||HXi+X2, 
where 

X, = -(divK®v^),vf), 
I, = (div(S(p^,D(v^))),vf). 
Using the Holder inequality and (13. lip , we deduce 

Xl < ||div(v^®V^)||2||vf||2<C||v^||3||Vv^||6||vf||2 

< ^llvfll^ + CllVv^lbllv^lhllV^v^ll^. 
o 

Next, using f lAip - (lA2p . we conclude 

X2 < II div (S(p^, D(v^))) II2 ||vf II2 < ^||vf 11^ + C (II Vp^ll^ + llv^liy . 
Finally, substituting the estimates for Xi and X2 into (15. 9p . we obtain 

(5.10) ||vf ||2 + C6j^\\v\''\\l < C(l + llVp^ll^ + ||Vv^||2||v^||2)(l + II VV^II^). 
Consequently, using (15. 8 p and the Gronwall inequality we conclude 



(5.11) sup ||VV^(t)||^+ / ||vf||^dt<C((5). 

iG{0,T) Jo 

Having (15.81) and (15. lip and using (12. 5p . we see that there is (v,p, S) such that 
(up to subsequence) 

^* V weakly* in L~(0,T; W^'^R^), 
vf ^ vt weakly in L'^{0, T; L'^{RY), 
^* p weakly* in L'^{0,T;W^^'^{R^)), 
S(j9^,D(v^)) -* S weakly* in L°^{0,T; L^{Ry^). 



^Notc that we use the fact that ||Vv^||2 = (Vv^, Vv^) = -(Av^, v^) < IjV^v 



'-^"l + lKlli- 
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It remains to pass to the limit in the equation f l5.5p . As usually, the most delicate 
terms to deal with are the nonlinear one. The classical tool is to use some com- 
pactness arguments in order to obtain strong convergence in the adequate function 
spaces. Hence, by using Aubin-Lions compactness lemma, we get (up to subse- 
quence) 

strongly in L\0,T;Wl;'^{Ry), 
a.e. in [0,T] x 
a.e. in [0,T] x R\ 

Consequently, using Proposition 14. we deduce 

-^p a.e. in [0,r] x R^. 

Using these convergence results it is easy to deduce that S = S{p, D(v)) and to set 
— )■ cxD in (15. 5 p and obtain a solution (v,p) to the equations (14. ip and (15. 4p . where 
(15. 4 p holds for all i/' with div?/' = 0. Finally using the Helmolhtz decomposition 
and the equation for the pressure (14. ip we finally obtain the validity of (15. 4p for all 
given ■?/'. The proof of Proposition 15.11 is complete. □ 



N 

V — )■ V 



Vv^ ^ Vv 



6. Limit 5 -> 

In this subsection we denote (y^jp^) a solution from Proposition 15.11 Our main 
goal in this section is to set 5 — )■ in the equation (15. ip and to get (12. 9p . 

First, using weak lower semicontinuity of norms and the Fatou lemma, it is easy 
to deduce from (15. 8p and from the estimate (13. 3p that 

sup ||v''(t)||2 + 2Ci /" [ f(^(|D(v^)|) + |v''|«dx) 

^-^ 0<t<T Jo Jrs V / 

+ 1 (5||VV||2 + C(5)||/||?,2) dt<C, 

for all q G [2, y). To obtain uniform estimate also on the pressure, we use Proposi- 
tion UT] to obtain the decomposition p^ = Pi + P2 such that for all q G (2, ^) and 
all s G [2, r'] there holds 

(6.2) r\\pi\\i + \\pi\\:<c{q,s). 

Jo 

Similarly, (15. lip implies that 

(6.3) sup II VV(t)||2 + r llvfll^ dt < C{5). 

t£{0,T) Jo 
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Therefore, using (16. ip . f l6.2p and the Aubin-Lions lemma, we can extract a not 
relabeled subsequence such that 





V 


V 


wpalrlv* in T,°°(C\ T- T?CU^'\^'\ 

Wccliviy 111 Ij lU, J. , iJ J I, 




V 


— ^ V 


weaKiy m L [0,i , vVi^^ } j, 


(6.6) 






weakly in Lif(0,T; W'^^'^^iW^f), 


(6.7) 




— V 


a.e. in (0,T) x 


(6.8) 


Pi 


^Pi 


weakly in L«(0, T; UiM?)) for all g G (1, — 

6 


(6.9) 


pi 


^P2 


weakly in L"(0, T; L"(M^)) for all s G [2,r']. 



These, estimates are sufficient to pass to the limit from (15. 4p to (12. 9 p if we show the 
point- wise convergence also for Vv^ and . 

To do it, we first derive some regularity estimates. For this purpose we observe 
first that v*^ is more regular. Indeed, solves 

(6.10) AV^ = F, 
where (after using (16. ip and (16. 3p ) 

F := div S{p\ D(v^)) - div (v^ ® v"^) - V/ - v\ G L'^{0, T; L'^{R^f). 

Thus, using (15. Sp . we get (for fixed 6) 

G L2(0,T; W^'^imy). 

Hence, we see that we can set := — Av*^ in (15. 4p . As the result, we obtain (after 
integration by parts and also using (lAip ) 

(6.11) + ^11 W(t)||^ + CM^') <Ji + J2, 

where 

Ir(v''): = j ^ (l + |D(v'^)|2)^ \B{Vv')\^dx, 

J2 : = (div (v^ ® v^) , Av^) . 
For ^2 we use the fact that div v"^ = to obtain 

j2 = -y2 f dkiy^'yd,ivydkivydx<c\\vv'\\i. 

Using the assumption (IA2p . we get 

Ji<lo [ (l+|D(v'^)|2)'^|D(Vv'^)||V/M^, 
<7o(l.(v^))^||V/||2. 
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Therefore, substituting the relations for J1-J2 into f l6.1ip . we get 
(6.12) 

^|l|Vv^(t)ll^+ 5\\W\\t)\\l+ CM^') < CWVv'Wl + 70 (I.(v^))'^' ||V/|b. 

To bound the right hand side of (16.121) . we use the estimate (14. 4p . Since, 70 < 
, we deduce that < Ci and therefore substituting (14. 4p into (I6.12p and 
using the Young and Korn inequahties, we finally deduce that (see [27] for details) 

(6.13) 1 |l|V^l2 + Clr(v^) < C (llD(v^)ll^ + llD(v^)ll^llVvi^) . 

Having the estimate (I6.13P at hand we proceed in a similar way as in [27| but now 
we need to be more careful since the behavior at infinity has to be discussed. 

6.1. Step 1. First,we recall two interpolation inequalities 

(6.14) Ml3<||-||^||-||L-" "=^^' 

(6.15) II ■ lis < II ■ 11(^11 ■ IlL"'^ ^ = ^- 
In what follows we also use the following notation 

|f I"* dx, Iklls > •= / kl"* 

'{\v\<l} '- J{\v\>l} 

The reason why we interpolate in (I6.14p -f l6.15p into L^"^ norm is the following lemma 
Lemma 6.1. The following estimate holds 

(6.16) I.(v^)>C||D(v^)ir3,,>. 
Proof. Since 

|V ((1 + |D(v'^)|2)3 - 1) |2 < C(l + |D(v^)|')^|D(Vv^)|2 

we obtain from standard embedding and from the simple inequality (1 + x)3 — 1 > 
Cxi, that is valid for some C > and all a; > 1, that 

ll|D(v^)|i|||> < C||(l + |D(v^)r)5 - 1||^ < CMv'). 

And (I6.16P immediately follows. □ 

To estimate the second term in the right hand side of (16.130 we compute 

l|D(v^)||^||Vvi^ < ||D(v^)||^,<||Vvi^ + ||D(v^)||^,>||Vvi^ 

< C(||D(v^)|||< + l)||Vvi^ + ||D(v^)0,>||Vv^||^ 



C(||D(v^)|||< + l)||Vvi^ + ||D(v^)||^;.;>«||D(v^)L^;^>||Vvi^ 

SEJei 2(1-/3) 
< C(||D(v^)|||< + l)||Vv^||^+(l.(v^)) ^ ||D(v^)||^|||Vv^||i 
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Applying the Young inequality, we conclude that 

i|D(v^)||^||Vvi^ < C(||D(v^)|||< + l)||Vvi^ + 5I.(v^) 

+ C||D(v^)||;;'^||Vv^|||^. 

2 

Finally, using the fact that ||D(v'')||r > < ||Vv^||2 we conclude 

l|D(v')ll3l|Vv1^ < C(||D(v^)||^^< + l)||Vvi^ + 5I.(v^) 

(6.17) 



+ C||D(v^)||:>||Vv^"^^- 



_4 

r 

12 



Next, for fixed i^i G (0, 1) we estimate the first term on the right hand side of 
( I6.13P similarly. Hence 

l|D(v^)||^ < CWVv'Wl + l|D(v^)||^;^>||D(v^)l|^;^^) 



< ||Vv'^||2+ (I,(v'^)) ^ ||D(v'^)||,'^(^-^)''||Vv^||^'^°. 

Finally, applying the Young inequality we observe that 

(6.18) ||D(v^)||^ < llVvi^ + 5I,(v^) + C||D(v^)||4||Vvif 

with A, B given as 

^^._3 (r-l)(-l + ^)(3r-2)r 

' -15r2 + 37r-27^r-18 + 18/i + 9//r2' 

(6.19) 

^ ~ 

-15r2 + 37r - 27 /ir - 18 + 18/i + 9/ir2 ■ 

Next, setting (note that for r G (|, 2), G (0, 1)) 

, . 3r2-llr + 6 
(6-20) , 

we observe that A = r and 

4(3 - r) 



(6.21) B 



3r — 5 



Since B > for all r G (1, 2) we obtain by combining fl6.17p . fl6.18p and fl6.13p 

that 

(6.22) |||Vvi^ + CI.(v^) < C{1 + ||D(v^)||;> + ||D(v^) |||<)(1 + HVv^^)^ 
with 

2(3 -r) 



(6.23) A := 



3r — 5 
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Finally, dividing fl6.22p by (1 + ||Vv^|||)'^, integrating w.r.t. time and using the 
estimate (16. ip . we finally deduce (see [29] for details) 

(6.24) [ I,(v^)(l + II Vv^||2)-^ dt<C. 

Jo 

6.2. Step 2. Here, we deduce how f l6.24p implies local compactness of D(v'^). 
Lemma 6.2. There is j ^ (0, 1) such that 

(6.25) [ (Ir(v'^))^ dt < C. 

Jo 

Proof. We can compute 

(I.(v^))^ = (I.(v^))^ (1 + II Vvi^)--^(1 + II Vvi^)-^ 

< C (I,(v^))^ (1 + ||Vvi^)-7A(i + ||D(v^)||2< + ||D(v^)|||>)^^ 

3r-2 3(2-r) 

and by using interpolation inequality || ■ II2 < || • ||r || ■ Hsr" and fl6.16p we deduce 
that 



(I,.(v^))^ < C {M^r')y (1 + ||Vvi^)-^^(l + ||D(v^)||^_< 

Sr^ 3(2-r) 

+ l|D(v^)IL,,l (i,.(v^)) ^ y\ 

Using the Young inequahty with coefficients ^ and ^3^, we continue as 
(I.(v^))^ < C^^^l^^ + C{1 + ||D(v^)|| -^) 

3r-2 -yX 3(2~r) j\ 

+ C||D(v^)||,^^ (I,(v^)) - . 

Finally, applying once again the Young inequality onto the last term and integrating 
w.r.t. time we see that we can choose 7 so small that by using f l6.24p and (16. ip we 
obtain Q 

The next lemma finally gives the desired compactness of the velocity gradient. 

Lemma 6.3. The following estimate holds 

2,J\r 



(6.26) / / iWrdx) dt<C{R). 

'0 \Jb{o,r) 



Proof. First, we can deduce that 

[ |D(vv^)r(ix 

Jb(o,fi) 



< I ((l + |D(v'^)n — |D(Vv^)n (l + |D(v'^)n 4 rfx 

B{0,R) 
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Using the Young inequality, we deduce 

/ \B{Vv^)\'dx< [ (l + |D(v'^)n'^|D(Vv^)|2(ix 

Jb(o,r) Jb(o,r) 



+ C I {l + \D{^')Y) dx. 

Jb(O.K) 



Finally, using 7-power, integrating w.r.t. time and using f l6.25p . (16. ip and the fact 
that we integrate over the domain with finite measure, we obtain fl6.26p . □ 

Finally, using the standard interpolation, (16. ip . f l6.26p and the Young inequality, 
we conclude 

iiv'^iT' < / iiv'^ir^^"'^^ iiv'^ir'^ < r 

\\l+a,r;B{0,R) - / 11^ II l,r;B(0,_R) II ^ H 2,r;B(0,_R) - 


provided that a ^ 1 is sufficiently small. Using the compact embedding W^^'^'^ into 
W^'^ we then conclude (up to subsequence) that 

D(v'^) ^ D(v) a.e. in (0, T) x W^. 

We skip the details here and we refer to [221 EZ] for a complete presentation. Con- 
sequently, we may use Proposition 14.11 and to finally deduce that 

/ a.e. in (0,r) x 



which finishes the proof of Theorem 12.21 
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